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Abstract 

A Bethe ansatz equation associated with the Lie superalgebra 
osp(l\2s) is studied. A thermodynamic Bethe ansatz (TBA) equation 
is derived by the string hypothesis. The high temperature limit of the 
entropy density is expressed in terms of the solution of the osp(l\2s) 
version of the Q-system. In particular for the fundamental represen- 
tation case, we also derive a TBA equation from the osp(l\2s) version 
of the T-system and the quantum transfer matrix method. This TBA 
equation is identical to the one from the string hypothesis. The cen- 
tral charge is expressed by the Rogers dilogarithmic function, and 
identified to s. 

1 Introduction 

Recently, thermodynamics of quantum integrable spin chains related to su- 
peralgebras received attentions. In particular, there are several papers fl], ||, 
[3], [|, [|, D, [7j on thermodynamic Bethe ansatz (TBA) equations related 
to sl(r\s). Namely, the TBA equation for the supersymmetric t — J model, 
which is related to the vector representation of s/(2|l), is proposed in Jl|. 
This is also generalized to the case of the fundamental representation of 
sZ(s| 1). The TBA equation for a supersymmetric extended Hubbard model, 
which is related to the fundamental representation of si(2|2), is proposed in 
|§. One can also find TBA equations beyond the fundamental representa- 
tions. The TBA equation of a generalized t — J model, which is based on a 
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tensor representation of s/(2|l), is proposed in |J. As is well known, sl(r\s) 
admits finite dimensional representations with continuous parameters. In 
fact the TBA equation of the model for electrons with generalized hopping 
terms and Hubbard interaction, which is related to one parameter family of 
four dimensional representations of s/(2|l), is proposed in Q. In the case 
of sl(r\s), the TBA equation with a continuous parameter is studied @ in 
relation with the continuum limit of integrable spin chains. These TBA 
equations are derived by the traditional string hypothesis R |10[. One can 



also derive the same TBA equations by the quantum transfer matrix (QTM) 
method |TJ, [T^, [□], [TJ], [T5|, |5j . In fact the TBA equations for the supersym- 
metric t — J model |J and the supersymmetric extended Hubbard model 
are also derived pj by the QTM method based on the sl(r\s) version of the 
T-system |16[ (a system of functional relations among transfer matrices of 
fusion models). 

On the other hand, study on TBA equations for osp(r\2s) case has be- 



gun quite recently in Refs. [LT], [18], in which we deal with the sim- 
plest osp(l\2) model from the point of view of the string hypothesis and 
the QTM method. The i?-matrix of the osp(r\2s) models are given as so- 
lutions of the graded Yang-Baxter equations |20|| , which have bosonic and 
fermionic degree of freedom. The i?-matrix of the osp(r\2s) models are pro- 
posed in |2~H| . Moreover they are constructed from the point of view of 
automorphisms of Lie superalgebras and classical graded Yang-Baxter equa- 



tions P^fl ; quantum supergroups [ 2"Bfl ; braid-monoid algebras [E^j. Some spe- 
cial examples are also examined from various context. See, for example, 
Refs. [^5[ |26], |27j. As for eigenvalue formulae of transfer matrices, osp{l\2) 
model and osp{l\2s) model are solved by the analytic Bethe ansatz; 
osp{2r — 1|2), osp{2\2s — 2), osp{2r — 2|2), osp{l\2s) models ^9f are solved by 
the algebraic Bethe ansatz. These works deal with models related to fun- 
damental representations. As for more complicated representations, we have 
executed an analytic Bethe ansatz |3{], [JTJ related to C(s) = osp{2\2s — 2) 
32] , B{r\s) = osp{2r + l|2s) and D{r\s) = osp{2r\2s) fl3~3| , in which 'su- 



per Yangian analogue of Young supertableaux' are proposed and T-systems 
among them are also found. The osp(r\2s) integrable spin chain is related 
to interesting physical problems, such as the loop model which is related to 
statistical properties of polymers and the fractional quantum Hall effect 
|35| , etc. So it is desirable to study the osp{r\2s) integrable spin chain beyond 
the osp{l\2) case. 

The purpose of this paper is to study TBA equations related to osp{l\2s) = 
B(0\s) by the string hypotheses [@, [TOf and the QTM method JTTJ, |TJ, |I§ [t§ 
pT5] , §|. In section 2, we assume a rather general Bethe ansatz equation (BAE), 
and derive the TBA equation by the string hypothesis. We expect this TBA 
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equation is the one for a higher spin analogue of the osp(l\2s) integrable 
spin chain in Ref. ||28|| . In fact, the high temperature limit of the free en- 



ergy is expressed by an appropriate solution of a difference equation called 
the Q-system. In particular for the fundamental representation case, we also 
derive the TBA equation from the QTM method in section 4. Namely, we 
transform the ID quantum system to the 2D classical system by the Suzuki- 



Trotter mapping [fLlJl , and define a QTM, which is a transfer matrix of an 
inhomogeneous vertex model. We consider the QTM in the context of the fu- 
sion hierarchy of the model, and derive the functional relations among 'fusion 
QTMs' from the osp(l\2s) version of the T-system|3B[]. After a dependant 
variable transformation, we obtain the F-system (which is expected to be a 
system of functional relations among a solution of the TBA equation) from 
the T-system. Finally we transform the F-system with certain analytical 
conditions (ANZC conditions: Analytic NonZero and Constant asymptotics 
for the large spectral parameter) into the TBA equation. Moreover we find 
that this TBA equation coincides with the one in section 2. This indicates 
the validity of the string hypothesis for the osp(l\2s) model. In section 5, we 
evaluate the low temperature asymptotics of the specific heat. We express 
the central charge by the Rogers dilogarithmic function and identify it as 
c = s. This coincides with the conjecture p8| by the root density method. 



2 String hypothesis 

For the last several decades, many people recognized that Bethe ansatz equa- 
tions (BAE) can be written in terms of the representation theoretical data of 
Lie algebras |36| or Lie superalgebras |2l], . So we assume, as our starting 
point, the following type of the osp(l\2s) Bethe ansatz equation on complex 



variables {v^}: 




(a) , _Lj>)\ s+1 n /,(<*), ,p 



J-j- Va(d)^fc Tj^J Q 1) 



i~A a) J £a l\Q,( d )(v ( k a) -IB ad ) 

where k G {1,2,..., M a }; a G {1, 2, ... , s}; a(d) = d for 1 < d < s; a(s + 
1) = s; t a = 1 for 1 < a < s - 1; t s = 2; B ad = 2S ad - 5 a4+1 - 
Q a (v) = rjj=i(w - ujfc ); M a e Z> ; M = N; N is the number of the lattice 
site, G C is an inhomogeneity parameter. In this section, we consider 
the case = 0. The parameter a expresses an effect of 'a peculiar two- 
body self-interaction for the root {f^}' |[29| , which originates from the odd 
simple root a s with (a s \a s ) ^ 0. The parameters (bf\ . . . , b^) in the 
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Figure 1: The Dynkin diagram for the Lie superalgebra B(0\s) = osp(l\2s) 
(s G Z>i): white circles denote even roots; a black circle denotes an odd root 
a s with (a s |a s ) 7^ 0. The Kac-Dynkin label b a of an irreducible representation 
of osp(l\2s) is depicted under each vertex. This irreducible representation 
is finite dimensional if and only if b a G Z>o for a G {1, 2, . . . , s — 1} and 
b s G 2Z> . 

BAE fl2.ip represent the Kac-Dynkin label [0 (see, Figure [I]). In this paper, 
we consider the case 

bf ] =b5 ap (l-5 ps )+2b5 ap 5 ps , (2.2) 

where b G Z>i, p G {1, 2, . . . , s}. Note that bj is always an even number. 
We write the irreducible representation of osp(l\2s) labeled by (|2.2|) as V^ v \ 



e a is a phase factor (|e | = 1). In particular, for p = b = 1 case, we have 128 ] 



(-1)N-M 2 jf a = 1 

(-l)M a -i-M a+1 jf aG {2,3,...,s-l} (2.3) 
(-1)^.-1-^. if a = s. 



Note that BAE (|2J) for p = 6 = 1 corresponds to the BAE [||, |9| of the 
2s + 1 state osp(l\2s) model [2l| |2^, E3| (the Hamiltonian of this model is 
given in Q4.2D ). 

Next we will present a function Tm\v) with a spectral parameter v G C, 
which is a candidate of an eigenvalue formula (DVF) of a transfer matrix 
for an osp(l\2s) vertex model of a fusion type for auxiliary space labeled 
by (a, m). We consider the case where the transfer matrix is defined as the 
ordinary trace of a monodromy matrix. The auxiliary space of the corre- 
sponding transfer matrix should be an irreducible finite dimensional module 
of the 'super Yangian Y(osp(l\2s)) '. As an osp(l\2s)-mo dule, is ex- 

pected to contain vjh* as 'top term': any other irreducible component in W$ 
has a highest weight lower than that of Vm ■ As for Wm \ the multiplicity of 
the irreducible components may be calculated by a Kirillov-Reshetikhin 



like formula in Ref. ||33|| . The quantum space that the transfer matrix acts on 
is the tensor product module (W 6 )® JV . We can derive an explicit expression 
of Tm ^ (v) by modifying the vacuum part of the DVF in Ref. EBJ so that the 
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vacuum part is compatible with the left hand side of the BAE (|2.1|): 



{d jk }j=lk=l 



-(m- a- 2j + 2k)), 



(2.4) 



where the summation is taken over d^ G {1, 2, . . . , s, 0, s, . . . , 2, 1} (1 -< 2 -< 
•■■-<s-<0-<s-<----<2-<l) such that djk z< dj+i^ and djk -< dj^+i- The 
functions {z(a;v)} are defined as 



z(a;v) = ip a {v) 

z(0;v) = i/j (v) 
z(a; v) = ih(v) 



Qa-i(v + %(a + l))Qa(v + %(a-2)) 
Q a -i(v + i(a-l))Q a (v + ia) 

for a E {1,2, ... , s}, 
Q a (v + i(s-l))Q a (v + i(s + 2)) 
Qs(v + ^(s + l))Q s (v + ^s) ' 
Qg^jv - i(a - 2s))Q a {v - |(g -2s- 3)) 
Q a -i(v - \ (a -2s- 2))Q a (v - §(o -2s- 1)) 
for a G {1,2,..., s}, 



(2.5) 



where Qo(' y ) := 1- The vacuum parts of {z(a; v)} are given as follows 



Ca 
Ca 

Ca 



for 1 ^ a ■< p, 



4>b{v+^ P -l))<t>b{v-^p-2s)) 
4> b {v+Up+l))<t> b {v-Up-2s-2)) 



for 



(2.6) 



for p ^ a ^ 1, 



where 0&(i>) = nlUi^ + 1 + 1 — 2&;)) 7V ; Ca is a phase factor (|£ | = 1): 
Ca = Ca = Ci^i^2 • • • e -i for a G {1,2,..., s}, Co = Ci^i^ In particular 

for p = b = 1 case, the phase factor Ca is given p8| as follows 



(-l)^ Ml if o = l 

(-l)^-x-w. if ae{2,3,...,s} 

1 if a = 

(-l)^-i-^ if aG {s,...,3,2} 



(2.7) 



Ca 



where a = a. The dress part of (|2.4|) is free of poles under the BAE ( |2.1|) . 
This is a requirement from the analytic Bethe ansatz pp| , |3"T[ |. Note that 
T^\v) for p = b = 1 corresponds to the DVF |2§ of the 2s + 1 state 
osp(l\2s) model. We may think of ( |2.4j) as osp(l\2s) version of the Bazhanov 
and Reshetikhin's eigenvalue formula p9f . 
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We shall consider the one dimensional counterpart of the above-mentioned 
fusion vertex model for the case (a,m) = (p,b), which is a higher spin ana- 
logue of the 2s + 1 state osp(l\2s) spin chain|2l|. The energy density of the 
corresponding system is defined as follows: 



— , (2.8) 



v= o N^( v M)2 + (kv 



k=l \ u k I T \2- 



where J is a real coupling constant. For simplicity, the function ( [2.6D is 
normalized so that the constant term of the energy density (|2.8|) vanishes. 



The extra signs (p.7|) do not appear in ( |2.8| ) explicitly. We adopt the following 
ordinary string solution for ( |2.1| ) i n the thermodynamic limit (cf. Ref. [[HJ): 

vS k + ^(rn + l-2j), (2.9) 

where m G Z>i; a G {1,2, ...,s}; j G {1,2, ...,m}; G {1,2, ... ,17^}; 

G R; nf is the number of color-a m-strings. Let pm(v) and (7m\v) be 
string and hole densities for each color a and length m. Substituting fl2.9| ) 
into ( |2.1|) , one can derive a relation between pm\v) and (7m (y) after some 
manipulation: 



mw(«) = + E E A ^V)> ( 2 - 10 ) 

1=1 d=l 



where the function ^ m ,b(v) is given by 

T /n i 3 f^ + i(m+l-2j + 6n 

27r ov [v + 7}(m + 1 — 2j — b) J 

and the operator is defined by 

A^l = 5 ad ([|/ - m|] + 2[|Z - m\ + 2] + ■ ■ ■ + 2[l +m - 2] + [I + m}) 
-I ad ([\l -m\ + l] + [\l-m\+3] + --- + [l + m- 1]), 

lad = 8a,d-l + $a,d+l + &ad&as- (2-12) 

Here the action of [m] on any function h(v) is defined as follows 
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where * is a convolution 

/oo 
dw^ mtl {v-w)h{w). (2.14) 
-oo 

In the thermodynamic limit, the energy density (|2.8| ) reduces to 

OO poo 

S = 2njJ2 dv* mjb (v)p%\v). (2.15) 

m=l 

The entropy density leads as follows 



S OO /»00 

s = k B E E / rfu i^ a) ^) + iog(pL a) (^) 

a=l m=l 



<7, 



(a), 



-pW («) log P W ( V ) - («) log o-W («) } , (2.16) 

where denotes the Boltzmann constant. Now we shall investigate the 
equilibrium state. From the condition — — = for the free energy density 

T = E — TS (T: temperature) and the relation ( J2.10 ), we obtain a TBA 
equation 

oo s 

log(l + y W («)) = 7*3 + E E A -1 + ( 2 - 17 ) 

Z=l d=l 

and the free energy density 

T =—aY.\ M^WMl+^W)- 1 ), (2.18) 

where Ym\v) := <Tm\v)/pm(v) and (3 := l/(ksT). We shall derive an 
alternative form of the TBA equation ( 2.17|) . Taking note on the relations 



oo 

ikv 



/oo 
dv^f m}b (v)e 
-oo 



• i / min(fe,m)A; \ / /t \i ; 
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we can rewrite fl2.17| ) as follows: 



1 via) l \ njpSapSmb 
log Y m {V) = _„t._.. + K * 



COsll 7Tf 



re =1 (i+(^ d) )- i ) / - 



(v), (2.20) 



where a G {1, 2, ... , s}, m G Z>i, Yq (u) := and the kernel is 

1 



K v) 



2 cosh TTV 



(2.21) 



Note that (|2.20|) reduces to the TBA equation in Refs. [|17|, [18|], if we set 
p — b — s — 1. 



3 High temperature limit 



We shall consider the high temperature limit of the free energy density ( [2.18 ) 
We assume that the function Ym\v ) is independent of v in the limit T —> oo 



In this case, the constant solution of ( p. 20 ) obeys the constant F-system 



where Yq := 0, a G {1, 2, . . . , s} and m G Z>i. The solution of this constant 
Y"-system is 



(3.2) 



where {Q^m} are the dependant variables of the Q-system: 

{Qt ] f = QtliQtli + fl(Q^) Iad i (3.3) 



<i=i 



where := 1, a G {1, 2, . . . , s} and m G Z>i. We can derive this Q-system 
by neglecting the vacuum part and the spectral parameter dependence of the 
T-system in Ref. [Q (see, ( |4.7|) ): a system of fusion relations among com- 
muting transfer matrices of the osp(l\2s) vertex model. The Q-system was 
proposed in Ref. |J8| as a system of difference equations among characters 



of modules of the Yangian Y(g) as g modules (g: simple Lie algebras). Our 

_____ /Q} ___ 

Q-system ( pj.3| ) is equivalent to the A K 2 J version of the Q-system . We 
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think that this coincidence originates from the correspondence between 
B^(0\s) and A^J ■ There are several similar correspondences between su- 
peralgebras and ordinary algebras. Thus, we expect that "different algebras 
with the same Q-system" also occurs to several other algebras. 

We can express the solution of the Q-system (|3.3| ) as a polynomial of the 
fundamental variables Qi, Qx, • • • , Q\ ■ Moreover, this polynomial has a 
determinant expression: 



where <2,- satisfies 



and a condition 



QL a) = det (Qa+i-j), (3.4) 



Qa = Q 2s -a+l (3.5) 



for a G Z <0 

Q a = { 1 for a = (3.6) 

Q[ a) for ae {l,2,...,s}. 

In our case, Q a is given as follows 



In this case, Q$ is the number of terms in DVF Q2.4|) , which is expected |33] 
to be the dimension of Wm\ In particular, coincides with the dimension 
of V^~\ For general case, Q^m is greater than or equal to the dimension of 



Vi a) . When solving ( [OOP , we assume Ym\±oo) coincide with the solution 



^2|) of the constant Y-system with (|3.4| )-([3T 



Substituting (|3~^ ) into (|2.18| ), we find that the high temperature limit 
of the entropy density is expressed in terms of a constant solution of the 
Q-system (|3.3| ): 



.7-. OO 

lim 5 = - lim - = k B V log(l + (yW)- 1 )™^ = k B \ogQ^. (3.8) 

T — >oo T — s-oo 1 ' * 

m=l 

In particular for p = b = 1 case, we have 

lim S = fcfllogQ^ = k B \og(2s + 1). (3.9) 

T — >oo 

This is the logarithm of the dimension of or This result is 

compatible with the number of the state per site of the osp(l\2s) integrable 
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spin chain. In closing this section, we note the fact that the solution ( |3. 2| ) of 
the constant y-system ( |3.1|) can also be written as 

Yi a) = -f f, 3.10 

a[g - a) 

where g — 2s + 1. By using ( [3.10] ), one can derive an explicit expression of 
a constant solution of the Q-system (|3.3|) 

n (g) = f (m + g)\m\ \ m yr \ (k + a)(k + g - a) \ k 

^ m \(m+a)\(m + g-a)\J W \ k(k + g) J ' l ' J 

which will provide a dimension formula of the module Wm^ ■ 



4 T-system and QTM method 



In this section, we introduce an integrable spin chain[T2^, [2Dj associated with 
the fundamental representation and derive the TBA equation from the 

point of view of the QTM method|Q], [E| p), 0, |j. A more detailed 
explanation of the QTM method in the case of osp(l\2) can be found in Ref. 



The i?-matrix[E2~ 



*M , |2P| of the model is given as 
2v 



R( v ) =I + vP 



2v-g 



-E, 



(4.1) 



where 



■l) p(a)pib) S ad S bc ; E cd 



ab 



a a b(a 1 )aj; a, b, c, d G {1, 2, . . . , s, 0, s, . . 
a is (2s + 1) x (2s + 1) anti-diagonal matrix whose non-zero elements are 
a ay a = 1 for a G {1, 2, . . . , s, 0} and a a ,a = —1 for a G {s, s — 1, . . . , 1}; 
a = a; p(a) = for a = 0; p(o) = 1 for a G {1, 2, . . . , s} U {s, . . . , 2, 1}. The 
Hamiltonian of the present model for periodic boundary condition is given 
by 



H 



k=l 



,k+l 

+ -E k 

,k+l 

9 



(4.2) 



where J is a coupling constant: J > and J < correspond to the ferro- 
magnetic and antiferromagnetic regimes respectively; L is the number of the 
lattice sites; Pk,k+i and Ek y k+i act nontrivially on the k th site and k + 1 th 
site. The QTM is defined as 



Tr.n^, 



a 2k A U + iv ) R a 2k -i,j( U - iv ) 



(4.3) 
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where R^(v) = R^(v); Rjk{v) = tk Rkj{v) {t^ is the transposition in the fc-th 
space); N is the Trotter number and assumed to even. By using the largest 
eigenvalue T x (un, 0) of the QTM (|4.3|) , the free energy density is expressed 
as 

T=-\ lim loglfV^O), (4.4) 

p N—*oo 

where = — ^M- . From now on, we abbreviate the parameter u in T^\u, v). 



One can obtain the eigenvalue formulae of the QTM ( |4.3|) by replacing the 
vacuum part (gl| of the DVF (fOj) for p = b = 1 with that of the QTM. 
One may interpret the QTM as a transfer matrix of an inhomogeneous vertex 
model. In our case, the inhomogeneity parameters in the BAE (|2.1| ) for 



p = b = 1 take the values: Wj = iud a i for j G 2Z>i; = (—iu + y)<5 a i 
for j G 2Z> + 1. The vacuum parts of (|2.4 ) are given as follows 

C ^^Tt?;^ for a=l, 
Mv) = { Ca<M?#» for 2 ± a d 2, (4.5) 
4>_{v+^i) tor a - X ' 



where 0±(f ) = (v ± m) 2 ; ( a is the one in ( |2.7| ). For a G {1, 2, . . . , s} and 
m G Z>i, we define a normalization function 

^ (v) = n 3 m = i m=i - ggf^)^ - g^±gp 

We reset T^\v)/Uin\v) to T^V), wh ere T^\v) is defined by Q. Since 
the dress part of the DVF Tm\v) is same as the row-to- row case, Tm\v) 
satisfies the following functional relation, which has essentially the same form 
as the osp(l\2s) T-system in Ref . [|33] . 

T^(V + 1 -)T%\V -'-)= T^ +1 (v)T^Uv) + T m~ 1 H V ) T m +1 \ V ) 

for a E 1, 2, . . . , s — 1, (4.7) 



where 



To a \v) = <f>-{v + - gi) for a G Z>i, 

T4°)(«) = ^(v-^Mv + j*) for mGZ> 1; (4.8) 
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For s = 1, gm(v)Tm\v) coincides with the function Tm\v) in Ref. fL8| . For 
m G Z>i, we define the F-functions: 

m ( J ^(^W 1}, 
yW(«) = n^^iT^n ■ (4-9) 

By using the T-system Q4.7| ), one can show that the F-functions satisfy the 
following y-system: 

y« {v + 1)^> ( „ _i) = (l+^W + ^M) , (4,0, 

2 2 ns=i(i+(^(«))- i ) i » d 

where Yq (t>) = 0? a ^ {1>2, and m G Z>i. A numerical analysis 
for finite iV, w, s indicates that a two-string solution (for every color) in the 
sector N = M x = M 2 = ■ ■ ■ = M s of the BAE (|2lf) provides the largest 
eigenvalue of the QTM (|4.3| ) at v = 0. Moreover, we expect the following 
conjecture is valid for this two-string solution (see Figure ||, Figure [| Figure 

§■ 

Conjecture 1 For small u (\u\ < lj and a G {1, 2, . . . , s} ; every zero o/ 
Tm\v) is located outside of the physical strip Imv G [— |, |]. 

Based on this conjecture, we shall establish the ANZC property in some 
domain for the F-functions (|4.9| ) to transform the Y~-system ( 4.10 ) to nonlin- 
ear integral equations. Here ANZC means Analytic NonZero and Constant 
asymptotics in the limit \v\ — > oo. One can show that the F-function has 
the following asymptotic value 

hl*W = ^, (4.1D 
|v|->t» a[g — a) 



which is identified to the solution (|3.10|) of the constant Y~-system ( |3.1|) . 
From the conjecture and (|4.11|) , we find that the functions 1 + Ym(v), 1 + 
(Y^iv))- 1 in the domain Imv G [-5, 8] (0 < 5 < 1) and Yt\v) for (a, m) 7^ 
(1, 1) in the domain Imr G [— |, |] (physical strip) have the ANZC property. 

On the other hand, Y^ l \v) has zeros of order N/2 at ±i(| — tt) if u > 
(J < 0), poles of order N/2 at ±«(| + m) if u < ( J > 0) in the physical 

strip. Then we must modify Y±(v) as 

± NS al S ml 

% o) 00 = Y^(v) jtanh |(„ + i(i ± «)) tanh - z(~ ± «))} ' (4.12) 
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1.4 



E 1.2 
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coloM ^ 
color2 □ 



-0.1 -0.08 -0.06 -0.04 -0.02 0.02 0.04 0.06 0.08 0.1 

Re 



Figure 2: Location of the roots of the BAE for osp(l|4) case (N = 6,u = 
0.05), which gives the largest eigenvalue of the QTM T^\v). Both color 1 
roots {'Wfc 1 ''} and color 2 roots {v^} form three two-strings. 
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I o 



m=1 
m=2 
m=3 



-0.1 



-0.08 -0.06 -0.04 



-0.02 





Re 



0.02 



0.04 



0.06 



0.08 



0.1 



Figure 3: Location of zeros of t£\v) for osp(l|4) case (m = 1, 2, 3, N 
6,-u = 0.05). The zeros recede from the physical strip \mv G [— \, \] as 
increases. 
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I o 



* x 



_l X L. 



m=1 + 
m=2 x 
m=3 * 



_l X L. 



-0.1 -0.08 -0.06 -0.04 -0.02 





Re 



0.02 0.04 0.06 0.08 0.1 



Figure 4: Location of zeros of Tm\v) for osp(l\A) case (m = 1,2,3, iV = 
6,-u = 0.05). The zeros recede from the physical strip Imv G [— |, |] as m 
increases. 
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tanh — (v + i) tanh — [v — i) = 1, 



where the sign ± is identical to that of — u. Taking note on the relation 

(4.13) 

one can modify the lhs of the ^-system ( [4.10[ ) as 

;i+^iW)(i+4 a i^)) 



2 n s d= i(i+(Yi d \v))-^ * 

for m G Z>! and a 6 {1,2, 



(4.14) 



s}. 



Now that the ANZC property has been established for the F-system, we can 
transform (|4.14|) into a system of nonlinear integral equations by a standard 
procedure. 



_N5 al 5 ml 



T- 



+K * log 



7T 



7T 



1 



log •jtanh —(v + i(— ± u)) tanh — (v — i(— ± u)) 



m=i(i+(^ 



-iy„ 



(4.15) 



where 3^ a ^(w) = 0, a G {1, 2, . . . , s} and m G Z>i. Substituting u = — A 
and taking the Trotter limit iV — > oo, we obtain the TBA equation ( |2.20|) for 
p = b = 1 . Taking note on the relations 

min(a,ci) 

Cad(^) = G k-rf|+2/-l( V ) ) 



/3J 



C ad (k) 



i=i 
4 



(2s+l-2a)7r 
4s+2 



2nv 
2s+l 



25 + 1 C o S ^t^^ +COSh;^ 



2s+l 



2s+l 



dwC a(i (v)e" 



-ikv 



^2C ac (k)D cd (k) = 5 ad} 

c=l 

Dcd(k) = 25 cd cosh - - J C(i , 

one can also rewrite this TBA equation as 
logY^(v) = 2-Kf3J5 ml G a (v) 



(4.16) 



+ y c aj * log ( 



(v), (4.17) 
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where Y 1 (v) = 0, a G {1,2, ... ,s} and m G Z>i. In contrast to ( |2.20|) , 
( [4.171 ) does not contain 1 + (Y^ (i*)) -1 which is not relevant to evaluate the 
central charge for the case J < (see, the next section). One can derive the 
following relation from ( [i.7|) for a = 1, (|4.9| ) and (|2.19|) . 

logT^^w) = log0_(f + i)(f)+(v — i) 

oo 

+ £*i,™*log(l + (Y«)- 1 )(<;). (4.18) 

m=l 

Taking the Trotter limit iV — > oo with u = — we find that JF = — i logT-L (0) 
with (|4.18|) coincides with the free energy density ( |2.18| ) for p = b = 1 (if we 
neglect the first term in the rhs of (|4.18|) coming from a different normaliza- 
tion of the DVF). One can also derive the following relation from ( |4.7| ) for 
m = 1, O and Q4TT6Q . 



loglfV) = log<Mi> + i)<Mu-i) +$^G a *log(l + r i (a)N 

a=l 



2e 2 sinh(A;tt) cos(kv) cosh(^-J-fc) 
£;cosh(^±±A;) 



+^ / dfc , ; Tirm ( 4 - 19 ) 



Taking the Trotter limit N — > oo with u = —4^, we obtain the free energy 
density without infinite sum. 

T = j{—?—(2\og2-^(—}—) + 4>(^^) 
\2s + lV ^ V 2s + 1 ; ^ v 2 + 4s ; 

5 POO 

-k B Tj2 / cbG» log(l + Y} a \v)), (4.20) 

a=l 

where ip(z) is the digamma function 

1>(z) = ^-\ogT(z). (4.21) 

dz 

The first term in the rhs of (|4.20|) for J = — 1 coincides with the grand state 



energy of the osp(l\2s) model in 28 



5 Central charge 



It was conjectured |28j that the osp(l\2s) model based on the fundamental 
representation is governed by the c = s conformal field theory. We shall 
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evaluate the central charge for the case p = b = 1 by using the result in the 
previous section and technique in |44, 45]. In this section, we assume J = — 1 
and ks = 1. We define a scaling function in the low temperature limit 

yS ± (v)=}imYi a \±(v + ^\og^)), (5.1) 

where 7 is a constant and Vp is the Fermi velocity [p8fl : Ujr = As ym\( v ) 

and behave in the same way, we set ym{v) = y^±(v). From the 

TBA equation Q4.17| ), ym\v) satisfies the following equation 

1 fa)/ \ 4w F 5 m i . v F a 

logjCW = e " F "sm - -) 

7 2 



115=1(1 + w 

where y^\v) := 0, a G {1,2, ...,s} and m G Z>i- The first term in rhs 
of ( |5.2|) is proportional to the a-th component r a of the Perron-Frobenius 
eigenvector of I ac : 

^ IqcT c = ( 2 cos - j r Q for a G {1, 2, . . . , s}, 

c=l ^ ^ ' 

r fl = Ain(^), £> a 2 = l. (5.3) 

' " " a=l 

We assume the DVF (f ) is expressed as a product of the ground state part 
T± 9s (v) and the finite temperature correction part T^^ n (v). In particular 
we can express the finite size correction part for small T by the function 

s 

logTl 1)fn (v) =J2 G a* log(l + Vf') W (5-4) 



a=l 



7T7 



s POO 

cosh(v F v)J2 sin (—) / dwe- ,,FU 'log(l + t/S a) («;)) + o(T). 

a=l 2 

Using the relation G a (— v) = G a (v), one can derive the following relation 
from 



A 2 ^ /*oo 

— E sin (T) / d^e— log(l + yS a) H) 



a=lm=l I \l+|/m ; (00)/ Vl+2/V(-Oo) 
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where L(x) is the Rogers dilogarithmic function 



L(x) 



dy 



logy log(l - y) 



i-y y 

For small T, the leading term of the specific heat 



for < x < 1. 



(5.6) 



C 



d_ 
df 



d 



T 2 --logT; 

dT s 1 



(0) 



(5.7) 



is proportional |46|, |47| to the central charge c 

TTCT 



C 



3vf 



+ o(T). 



Then we obtain 
6 



7T- 



EE £ 



a=l m=l 



1 + ^(00) 



1 + ym (-00) 



(5.J 



(5.9) 



This expression is independent of the choices of the constant 7 and widely 
seen in the model related to rank s algebras |48|, |45|. To proceed further, we 



have to evaluate the limit ym (±00). At first, we expect lim^^oo?/: 
given by (|4.11| ). On the other hand, the divergence from the first term in rhs 
of (|5.2|) in the limit v — > —00 is expected to be canceled by lhs if 



(«), 



V) IS 



-00. 



(5.10) 



Then ( |5.2|) in the limit v — > —00 will be valid if ym\— 00) for m G Z>2 
obey the constant F-system ( p.l|) . Thus we expect ym\— 00) is given as the 
solution ( |3. 10| ) of ( |3.1| ) with a shift m — > m — 1: 



lim ?/ 



(a)/ 



(m — l)(g + to — 1) 
a(g - a) 



(5.11) 



Using the relations (|OD , (|4TT| ), ( gH|) and L(l) = ^, L(0) = 0, we find 
c = s. This result is consistent with the conjecture by the root density 
method. 



6 Discussion 



In this paper, we have proposed the TBA equation for osp(l\2s). It will 
be an interesting future problem to study the TBA equation based on a 
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more general orthosymplectic Lie superalgebra osp(r\2s) (r > 1). As for the 
string hypothesis, we will need to consider more complicated strings than 
( |2.9| ) as sl(r\s) case[|I], Q [7J. On the other hand, we have only the T- 
system for tensor-like representations fl33| . To construct complete set of the 



T-system which is relevant for the QTM method, we have to treat spinorial 
representations. 



QTMs which act on Z 2 graded vector space are proposed in [^. |5"0 |. 
In particular in [fjOfl , two types of QTMs, which give the same partition 
function in the Trotter limit iV — > oo are proposed. The one is defined as 
the supertrace of a monodromy matrix of an inhomogeneous graded vertex 
model. In this case, the sign of the BAE (|2.1| ) will have the form e a = 
(— l) deg ( a<j ) ( deg(a a ) = for the even roots; deg(a s ) = 1 for the odd root). 
The other is defined as the ordinary trace of the monodromy matrix. In this 
case, the sign of the BAE ( |2.1|) will disappear (e a = 1). As far as the ground 
state case is concerned, the difference among these formulation of the QTM 
and the ones in this paper may not be important. 

A fermionic formula for A^J coming from a different labeling of the 



Dynkin diagram is proposed in |pl| . Taking account of the correspondence 
j43fl between A^} and B^(0\s), this formula may also be relevant to osp(l\2s) 



case. 
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Note Added 



After the submission of this paper, a preprint |52] appears, where TBA equa- 
tions related to osp(l\2s) models are discussed from the point of view of 
integrable field theories. 
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